Global Neutrino Data Analysis and the Quest to Pin Down 
sin ^13 in Different Mixing Matrix Parametrizations 
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Hints for sizable sin 613 have been reported in earlier global neutrino oscillation data analyses as 
well as will be reported in this work, and quite recently by the Double Chooz experiment. However, 
as we enter the era of precision neutrino oscillation experiments, terms linear in sin ^13 will no longer 
be negligible, and its sign would affect the extraction of other oscillation parameters. The sign of 
sin 613 also plays a crucial role in the determination of the CP- violating phase. In this work we show 
that by adopting an alternative parametrization for the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) 
mixing matrix, one already has a chance to infer the sign of each mixing angle in the conventional 
parametrization using existing global neutrino data. A weak preference for negative sin 9is is found. 
In particular, the solar data suggest that sin 613 > while all other data the opposite. This leads 
to the speculation on whether the Mikheyev-Smirnov-Wolfenstein (MSW) effect is responsible. In 
this work we found that in the new mixing matrix parametrization, the 68% CL constraints on 
the three mixing angles are comparable to those estimated in the conventional parametrization 
adopted in the literature. Owing to the strong correlations among the three mixing angles in the 
new parametrization, the advantages of doing the global neutrino oscillation analysis using data 
from past, current, and near future neutrino oscillation experiments shall become manifest. 

PACS numbers: 14.60. Pq, 12.15.Ff, ll.30.Er 
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I. INTRODUCTION 



The experiments involving solar, reactor, atmospheric 
and accelerator neutrinos have established a picture of 
neutrino oscillations caused by non-zero neutrino masses 
and mixing among different neutrino flavors (see e.g. 
Ref. [l|). The phenomenology of neutrino oscillations 
arising from the mismatch between the weak and the 
mass eigenstates can be described by the Pontecorvo- 
Maki-Nakagawa-Sakata (PMNS) mixing matrix d Q. 
This matrix can be parametrized in various ways as seen 
in the literature (e.g., 043)- 
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All experimental data except those from the LSND [l3| 
can be well described assuming three active neutrinos. In 
the case of Dirac neutrinos, the 3x3 unitary mixing ma- 
trix is characterized by three Euler angles and one phys- 
ical phase, and can be expressed as a product of three 
rotation matrices. The physical phase can be responsible 
for CP violation in the neutrino sector. As mentioned in 
Ref. the CP-violating phase can be associated with 
the sine or cosine of any mixing angle or with the iden- 
tity entry in any of the three rotation matrices. In this 
work, we will follow the standard Cabibbo-Kobayashi- 
Maskawa (CKM) matrix [ll| for the assignment of the 
CP-violating phase. 
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where 9ij and Sep are the mixing angles and the CP- 
violating phase, respectively, = cos 9ij and Sij = 
sinOij. As the combination R23W13R12 is the stan- 
dard choice for describing the quark mixing, it has been 
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adopted to be the conventional parametrization for the 
mixing matrix for Dirac neutrinos. This choice of the 
parametrization was actually made prior to the era when 
neutrino oscillation data became available. Later, the so- 
lar neutrino experiments fl^-flBI . reactor [igI . [Tt} , long- 
baseline (LBL) accelerator |18l - [20ll. as well as the at- 
mospheric neutrino experiments |2ll423l | all choose this 
parametrization to present their results. Using exis ting 
global neutrino data, most phenomenology works |24j3l| 
also employ this parametrization to determine neutrino 
oscillation parameters. 

In the conventional parametrization, the three mix- 
ing angles happen to nearly decouple for solar, atmo- 
spheric/accelerator, and reactor neutrino oscillation ex- 
periments due to the different neutrino energies and 
traveling distances involved. The 9i2 and 623 angles 
are well determined by solar and atmospheric experi- 
ments, respectively. Current-generation short-baseline 
reactor exp eriments such as the Daya Bay [32j . Double 
Chooz j33j |. RENO [S^ and the Angra [s^l experiments 
are exploited to pin down the yet unknown 6*13 value. 
Non-zero or sizable 613 values are pre dicted by many 
neutrino mass models (see e.g. Ref. [3g| for a nice com- 
pilation of existing model predictions) , and supported by 
global neutrino data analyses (see e.g. Ref |24| - [26l |29| | 
and this work). Very recently, the Double Chooz ex- 
periment 37] has reported their preliminary results of 
sin^ 2013 = 0.085 ± 0.051 (68% CL). When this is con- 
firmed in the future by other experiments with even bet- 
ter sensitivities, it will be good news for near- future ex- 
periments such as T2K [H, NOi^A T2HK [40|, 
T2KK ^ij. Their goal of measuring the CP-violating 
effect will be more reachable. 

This is not the end of the story for the neutrino os- 
cillation community. Another issue is the determination 
of the sign of sin^ia. As we enter the era of precision 
neutrino oscillation experiments, terms linear in sin^ia 
may no longer be negligible in fitting the mixing angles. 
From the Jarlskog invariant quantity [13, of CP vio- 
lation, one sees that the sign of sin 6*13 would also have 
an impact on the Sep determination. However, in the 
conventional parametrization, the reactor and solar ex- 
periments are only sensitive to sin^ 6*13 . For LBL accel- 
erator or atmospheric experiments, the observable has 
terms linear in sin ^13, but only sub-dominant. Despite 
the unfavorable situation one faces when working in the 
conventional parametrization, a first attempt to deter- 
mine the sign of sin 6*13 has been made in Ref. [131 using 
the LBL accelerator, CHOOZ, and atmospheric neutrino 
data. 

In fact the conventional parametrization is not the 
only way to establish the mixing matrix. As proposed 
by several authors [343 , the corresponding mixing pa- 
rameters may be more accessible, without sacrificing ac- 
curacy, in other parametrizations. In this work we fol- 
low the approach of Ref. [9] to explore this possibility. 
Besides the conventional parametrization, we perform a 
global neutrino oscillation analysis adopting an alterna- 



tive parametrization, R13W12R23, in which the observ- 
ables have leading terms linear in any of the three mix- 
ing angles. Equipped with the analysis results obtained 
in these two parametrizations, we will be eligible to ad- 
dress a couple of issues. Can one determine the sign of 
sin0i3? Are there other parametrizations which can pro- 
vide comparable sensitivities in extracting neutrino oscil- 
lation parameters as the conventional one? How are the 
three mixing angles correlated with each other therein? 
Do the matter (MSW) effects [H, SI] affect the sign of 
sin 6'i3? 

This paper is organized as follows. In Section |lT] we 
briefly describe our data fitting procedure in the conven- 
tional parametrization (to be denoted by A). In Sec- 
tion UTT] we present and discuss our results obtained from 
a similar analysis done in an alternative parametrization 
(to be denoted by D). Section ITVl gives our summary and 
outlook. Individual analysis approaches for each neu- 
trino oscillation experiment are detailed in Appendix 1X1 
through |D] 



II. ANALYSIS 

Table I in Ref. [^ indicates that the predicted neu- 
trino mixing angles in the parametrizations R23W13R12, 
R23W12R13, and R13W23R12 come up with similar 
values, while those in R13W12R23, R12W23R13, and 
-Ri2W^i3-R23 end up with similar values. We thus preform 
a global fit of the neutrino mixing parameters using pa- 
rameterizations R23W13R12 and R13W12R23, taking ad- 
vantage of the fact that these two parameterizations will 
have dissimilar outcomes. Following Ref. [^, we denote 
them as the A and V 'representations' respectively, where 
representation A corresponds to the conventional mix- 
ing matrix parametrization. In what follows, notations 
{A-^mf^, e{^) and (A'^mfj, Of^) stand for the oscillation 
parameters directly determined using representations A 
and V, respectively. On the other hand, the notation 6'^'^ 
symbolizes the three mixing angles that are initially ex- 
tracted from representation T) and then translated to rep- 
resentation A by the transformation method described in 
Ref. . The transformation from representation V to A 
is outlined in Appendix [E] 

The global neutrino data used in this work include 
(i) solar data from rates measured in chlorine [i^ and 
gallium [13] experiments, the rate of ^Be solar neutri- 
nos measured in the Borexino (48j experiment, Super- 
Kamiokande (SK) phase I & II day/night spectra Il4 1 1 
and SNO phases I & II j/g survival probability ll 
reactor data from KamLAND [Tiaiid CHOOZ [l 
LBL accelerator data from K2k|i1] and MINOS dis- 
appearance channel ^] and Ve appearance channel [20j : 
and (iv) atmospheric data from SK phase I [23j and 
SNO [11. 

In this work, the re-evaluated Df, flux from nuclear 
power plants [43| is not taken into consideration for the 
reactor data. In addition, the large number of bins in the 
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atmospheric data of SK phases II and III, as weh as the 
lack of information, prevents us from reproducing their 
results. Therefore we do not include the SK-II and SK-III 
atmospheric data in our work. Since our purpose is to in- 
vestigate the neutrino mixing phenomenology in different 
parametrizations using the same data sets and analysis 
conditions, the absence of these factors should not have 
any impact on the conclusions of this work. Described 
below is our analysis of the global neutrino data. 

Appendix |X] outlines the analyses for each solar ex- 
periment employed in this work which includes chlorine, 
gallium, Borexino, SK and SNO. We use the Bahcall so- 
lar neutrino spectra [50| except for the ®B neutrino spec- 
trum, which is from Ref. |5l| . Neutrino survival prob- 
abilities in the Sun are estimated using the BS05(OP) 
model [s^] for the neutrino production rates at different 
solar radii. We do not use the more recent BPS09(GS) or 
BPS09(AGSS09) models H because of the conservative 
model uncertainties. For estimating neutrino survival 
probabilities inside the Earth, the Earth density profile 
of PEM-C Q (rather than PREM [H) is adopted in 
this work. It assumes the continental crust for the outer 
most layer of the Earth where the solar neutrinos enter 
the detectors. In our analysis, the flux of ®B solar neutri- 
nos is a nuisance parameter while the fluxes of other solar 
neutrinos are taken from the BS05 model [s^]. Neutrino 
oscillation probabilities are calculated using the adiabatic 
approximation [56]. It has been verified to yield equiva- 
lent results as the numerical calculation does for several 
sets of neutrino oscillation parameters in the LMA re- 
gion. 

Our analyses of oscillation parameter fitting using the 
KamLAND and the CHOOZ oscillation events are de- 
lineated in Appendix |B] Owing to the short distances 
between the source and the detector, matter effects are 
generally insignificant and oscillation probabilities in vac- 
uum suffice. This is also the case for accelerator neutrino 
analyses. Appendix [C] describes the analyses of oscilla- 
tion parameter fitting using the K2K and the MINOS v^J^ 
disappearance channel as well as the MINOS i^e appear- 
ance channel. 

Atmospheric neutrino data from the SK [2^ and the 
SNO p^ l experiments are included in our analysis as de- 
scribed in Appendix |DJ We first employ the NUANCE 
package [57| to simulate atmospheric neutrino events as- 
suming no oscillation effects. Neutrino oscillations in the 
atmosphere and inside the Earth are then incorporated 
by the "weighting factors", Eq. (jDl|) . with the matter ef- 
fects taken into account by following the prescription of 
Ref. [58| . We apply similar criteria and cuts on the kine- 
matics of the simulated events so that we achieve the 
same selection efficiency as Ref. [1^ does for SK atmo- 
spheric data and as Ref. [12] does for SNO atmospheric 
data. 

We first perform a global analysis using representa- 
tion A in order to compare with existing results from the 
literature. As the current global neutrino data is insen- 
sitive to (Sep, we have selected 5cp = 0. According to the 



study presented in Ref. [!|, this choice also turns out to 
be adequate for representation T> since a rephasing still 
gives (5cp = therein. Table U summarizes our best-fit 
results for the oscillation parameters from various three- 
flavor analyses in this case. Assuming CPT invariance, 
the results from the KamLAND reactor experiment and 
all solar experiments are combined in the function: 
xLi+KL = xLiar + XkL' Wc minimize X^^^i+kl ^i*!^ re- 
spect to the oscillation parameters A'^m2i, tan^ 6*5^, and 
sin^ 6*5^, as well as the ^B solar neutrino flux in a three- 
flavor hypothesis. The other two parameters are fixed at 
A-^m|2 = 2.4 X 10-3 eV^ and sin^ = 0.5. The best- 
fit values are found to agree with the results presented 
in Refs. [H, [l^. Figure [T] shows the allowed regions in 
the planes of (tan^^^^, A-^m^i), (sin^Sjjj, A-^TO2i), and 
(tan^Sj^, sin^ Bf^. 

The x^ functions for K2K, MINOS t'^ disappearance, 
MINOS i/e appearance, CHOOZ, as well as SK and SNO 
atmospheric neutrino data are also summed together: 
xL+CHOOZ+atn. - xLc+X^HOOZ+xL„>. By minimizing 
this x^ with respect to the oscillation parameters A-^to§2, 
sin^ei^g, and sin^ 6*5^3 while fixing A-^to^i = 7.67 x 10"^ 
eV^ and tan^ Bf^ = 0.427, and assuming a normal mass 
hierarchy, the best-fit values are found to be consistent 
with the most recent results in Ref. [l^, . Figure [5] 
shows the likelihood contours in the planes of (sin^ 6*23, 
A-^m22), (sin^ e^^ , A'4m22), and {^m" B^^, sin^S^^). 

Finally, for individual data sample mentioned above, 
i.e. the solar-only, KamLAND, solar -|-KamL AND, 
and accelerator -|-CHOOZ+atmospheric, the differences 
(Ax^) between x^ and the minimum Xmin ^ ^ function of 
sin^ 6^^ are obtained by marginalizing the other two os- 
cillation parameters. Our results are shown in Fig. [3] for 

^Xsolar ' ^XkL' ^XsoH-KL AXacc-l-CHOOZ-l-atm- The 

global result is then 

^Xglobal(s"^^ ^13) = ^Xsol+KL + ^Xacc+CHOOZ+atm • 

(3) 

Our result for sin is consistent with those from other 
global neutrino data analyses (e.g., El, [ll H-il] ) • We 
also find a weak hint for non-zero sin B^^ at 95% CL, but 
it is not easy to determine the sign of sin^^^. 



III. RESULTS AND DISCUSSIONS 

Changing to representation V and fixing S^p = 0, we 
perform a similar x^ analysis as that in previous sec- 
tion using exactly the same data sets. After determin- 
ing the oscillation parameters A^nifj and Bj^, with the 
transformation strategy described in AppendixlEl we also 
perform a translation of the best-fit mixing angles back 
to representation A, denoted by Bj^-^. Both results are 
given in Table [TTl for a comparison with the corresponding 
A-^mfj and B^ retrieved directly in representation A. 

We found that in representation P, indi- 
vidual data sets, i.e. the global solar data, 
KamLAND, solar-|-KamLAND, and the combined 
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TABLE 1: Summary of our best-fit oscillation parameters, A'^m^j and 6^, in representation A, where the mxing angles, 0^, are 
in unit of degrees. In our analysis we only fitted sin^ ^23 in the range below 0.5. In order to compare the 68% CL constraints 
determined in representations A and T>, we assume the upper and lower bounds from accelerator+CHOOZ+atmospheric data 
are about the same. 
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TABLE II: Summary of our best-fit oscillation parameters, A^niij and 9j^, in representation ©, and the translated mixing 
angles, 9f^^ in representation A, where all of the mixing angles are in unit of degrees. Note that the 68% CL constraints of 
the best-fit results from the second minimum is taken at Ax' = x' ~ Xmin = 2.3 with Xmin being the x' value at the second 
minimum. The uncertainties of 6'^"^ are listed at the last row. 



accelerator+CHOOZ+atmospheric data, cannot con- 
strain the three mixing angles well. However, when 
combined altogether, they infer bounds at 68% confi- 
dence level (CL) on the three mixing angles Ofj which 
are comparable to those on 0^ (cf. Table HIl) . except for 
9f2- We therefore give only the bounds obtained from 
the combined global neutrino data. 

Using the global solar data, we find that the best-fit 
values of OV^-^ and 9V^ are consistent with those of 9-A 



and 6*]^, respectively. In representation A, the solar data 
are not sensitive to ^23 1 thus one should not naively com- 
pare the translated best-fit result of 0^^^ ^ 50° to the 
input value of 6*^3 = 45°. The KamLAND data are insen- 
sitive to 0^ too and have less constraining power on ■ 
These two facts explain why the translated values of 9^!/^ 
and 6*13^ are very different from the expected values, 9^ 
and 923- Surprisingly, the translated best-fit 6*^-^ turns 
out to be almost the same as 9^-^ in magnitude, but with 
an opposite sign. Like the KamLAND-only results, the 
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FIG. 1: Results of the solar and KamLAND oscillation pa- 
rameter analysis in a three-flavor hypothesis performed in 
representation A, where A'^m32 
0.5 are assumed. 



FIG. 2: Three-flavor oscillation parameter results in repre- 
sentation A for the combined analysis of the LBL accelerator, 
CHOOZ, and atmospheric data sets, where the parameters 
A-^mii = 7.67 X 10"^ eV^ and tan^ 6*5^ = 0.427 are assumed. 



combined solar and KamLAND results produce best-fit 
values where 9f/ - 6^^ and \9f^^\ ~ 1 6*5^3 1, and Of^-^ < 0. 

Similarly, combining the LBL accelerator, CHOOZ, 
and atmospheric neutrino data, the translated best- 
fit and \9f^\ are consistent with the best- 
fit and 16*5^1, respectively. The combined 
accelerator -|-CHOOZ-|-atmospheric data are not sensitive 
to 6-^, so a comparison of the best-fit value of 9^2^ ^ 45° 



with the input value of 6-^ = 33.2° is not necessary. In- 
terestingly for this case, it was found that the best fit 
9f^ < as well. This finding has been discussed in 
Ref. [23] using a similar data sample. 

In spite of the poor constraints on 9j^, the parameter 
A^m2i determined by analyzing the solar, KamLAND 
and solar -l-KamL AND data, as well as 77132 from the 
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FIG. 3: Ax'^ as a function of sin^ for various combina- 
tions of oscillation data where the undisplayed parameters 
are floated. 

combined analysis of the LBL accelerator, CHOOZ and 
the atmospheric data agree well with those directly ex- 
tracted in representation A. Equipped with the above 
individual calculation, we perform a global analysis 
straightforwardly via 

^Xglobal(^12i ^23' ^13 ) = ^Xsol+KL + acc+CHOOZ+atin ' 

(4) 

Here we fix A'^m^i = 7.6 x 10"^ eV^ and A^m§2 = 
2.4 X 10"'^ eV^, the best-fit values determined from the 
solar -|-KamLAND and the acc-|-CHOOZ-f atm analysis, 
respectively. Our results at 68% CL are given in Table HIl 
and Fig. 2] shows the projected likelihood contour plots 
in the planes of (^fa^^Fs), (^is^^Fs), and (^fa^^Fs)- As 
already advertised in the beginning of this section, we 
find that the 68% CL constraints on 6*^ are comparable 
to those on 6*^, with the exception of 



Sign of sin6'5lj 

As we enter the era of precision neutrino oscillation 
experiments, terms linear in sin^ja can no longer be ne- 
glected. As pointed out in Ref. [23|, it will be necessary 
to perform a parameter fitting in the sin^jlj < regime 
as well. The sign of sin^?i3 also plays the decisive role 
to the 8cp determination. This can be seen in the Jarl- 
skog invariant quantity of CP violation ^] . For any 
representation 7^, it is defined as 

J = sin 29^^ sin 261^3 sin 261^3 cos 61,^ sin j;?^ , (5) 

where 0^ is the mixing angle seated in the middle of 
the three rotation matrices in the mixing matrix. For 
representation H = ^ or 2?, both sin 20^ and sin 26*^ are 
positive. However, for the existing global neutrino data, 
the leading terms of the observables in representation A 
are linear in sin^ ^13 while the terms linear in sin ^13 are in 



the sub-leading terms. This makes difficult to determine 
the sign of sin ^13 in representation A. 

An interesting feature is seen when analyzing neutrino 
data in representation T). Two local minima exist in the 
global x^, one corresponding to 6'^'^ < (the global min- 
imum) and the other one d^-^ > (the second minimum). 
Apart from this, both minima have very similar values 
of 9f2^, ^23"^, and l^is'^l, which also agree well with the 
best-fit 9^2, ^23 > and \ef^\ (cf. Table IJ), respectively. This 
suggests that by doing global neutrino data analysis in 
representation 2?, one has a chance to identify the sign of 
sin 0^3^. In the current case, the x^ values for the global 
minimum and the second minimum differ by ~2. While 
this preference for negative sin0i3 is weak, the inclusion 
of data from the current and the near-future neutrino 
experiments may help strengthen the evidence. 

Compared among the various data samples used, the 
global solar data point to 9f^ > while all the other 
data 6*^3^ < 0. We do not have the explanation for this 
difference yet. One may investigate whether the matter 
(MSW) effects gl IH influence the sign of 9f/' during 
the propagation of solar neutrinos through the Sun and 
the Earth. Atmospheric neutrinos propagating through 
the Earth are also subject to the MSW effects (see Ap- 
pendix |D]), but to a much less degree. This issue will be 
further studied in a separate work. 



Error Correlations 

Another issue to address is how any two mixing angles 
correlate with each other. In principle, the 68% CL con- 
straints on 9j^-^ should be expected to be the same as 
those on 9^ since they are extracted from the same data 
sample. By virtue of this fact, one can estimate the corre- 
lation coefficients between any two mixing angles in each 
representation. For two different representations, say Q 
and v., any mixing angle in Q, 0^^,, can be expressed as 
a function of the three mixing angles, 9jj, in H as: 

^ab = ^(^12i^'23'^'l3) ■ 

By applying the error propagation 

the correlation coefficients p can be analytically solved 
using Eq. ([5]). With the best-fit results from the global 
minimum in 2? and the corresponding translated results 
in A, along with the average of upper and lower bounds 
of uncertainties from 9Yj and 6*^^, one can estimate the 
correlation coefficients p of any two mixing angles in rep- 
resentations A and P, respectively. Two situations are 
compared, (i) the 68% CL constraints taken at Ax^ = 2.3 
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68.27% C.L 

95.45% C.L 

99.73% C.L. 

X best-fit 



Rep. V 




9,2 (degree; 



Rep. p{6^2,92z) p{623,ei3) p{e^2,eri) 
Case (i): Ax^ = 2.3 



for all Mfj and AOfj 



A 
V 



0.288 
-0.854 



0.086 
-0.114 



0.290 
0.623 
-X — 



Case (ii): Ax =1.0 for A6I13 
but Ax^ = 2.3 for the rest 



A 
V 



1.861 
-1.082 



-0.554 
-0.409 



3.000 
0.632 



TABLE III: The estimated correlation coefficients of any two 
mixing angles in representations A and T), where the 68% 
CL constraints on and 613 are assumed by Ax'^ in two 
situations. 



Rep. V 


68.27% C.L 

95.45% C.L 

99.73% C.L 

X best-fit 
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FIG. 4: Three-flavor hypothesis mixing angle results for the 
combined analysis of global neutrino data in representation 
T). The neutrino mass differences are fixed to the values 
A^m^i = 7.6 X 10"^ eV^ and A'^miz = 2.4 x 10"^ eV^ 
and the normal mass hierarchy is assumed. 



for all (or 6lg-^) and Of^; and (li) the 68% CL con- 
straints taken at Ax^ = 1.0 for 6']^ while keeping the 
rest at Ax^ = 2.3. The results of the estimated correla- 
tion coefficients are presented in Table IIIII 

It is interesting to note that the correlations between 
any two mixing angles in representations A and V are 



bigger if case (ii) is applied. The outcome of the correla- 
tions for case (ii) seems to conflict with existing neutrino 
data which indicate the three mixing angles in represen- 
tation A are nearly decoupled. However, the results of 
the correlations in case (i) are more consistent with the 
nearly-decoupled feature among the three mixing angles 
in representation A. This suggests that the 68% CL con- 
straints on seem to require Ax^ = 2.3 rather than 
Ax^ = 1.0 since are also correlated with the other 
mixing angles (cf. Table HlH) . In addition, the correla- 
tions among the three mixing angles in representation A 
are not completely zero though small, indicating that the 
three mixing angles in representation A are not really de- 
coupled. Nevertheless, the correlations among the three 
mixing angles in representation A are found to be smaller 
than those in T), as expected. With the estimated p^, 
the 68% CL constraints on 6'^-^ for tlie second minimum 
are presented in Table [III 



IV. SUMMARY AND OUTLOOK 

We have performed global neutrino oscillation data 
analyses in two representations for the neutrino mix- 
ing matrix. Individual data samples we used include 
those from solar, KamLAND, long-baseline accelerator, 
CHOOZ and atmospheric experiments. We found that 
individually they do not constrain the three mixing an- 
gles in representation V as well as those in representation 
A. However, the 68% CL constraints on fff^ retrieved 
from combined global neutrino data in representation T) 
are almost as good as those on 6*^ determined in repre- 
sentation A, with the exception of OY2 which is a little bit 
worse than By^ by ^1.5°. Nevertheless, the best-fit Ofj re- 
sults, as expected, turn out to be significantly large. This 
results provide a higher sensitivity to the CP-violating 
phase determination as can be seen from the Jarlskog 
invariant quantity. In addition, the translated best-fit 
result of each 9f!j^ from representation I? to ^ are found 
to be in very good agreement with the corresponding 
angle, 0^, directly retrieved from representation A. In 
other words, 6*^-^, 9^^, and |6'J^-^| are respectively con- 
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sistent with determined from the combined resuh of 
solar+KamLAND, ^33 obtained from the combined re- 
sult of the accelerator+CHOOZ+atmospheric data sets, 
and 1^5^ I extracted from the combined resuh of global 
neutrino data in representation A. 

As we enter the era of precision neutrino oscillation ex- 
periments, observables with terms linear in sin^ia may 
be no longer negligible in fitting the mixing angles. The 
sign of sin plays a key role in the determination of 
the sign of 5cp, and both sin^ia and Sep are important in 
the establishment of the mixing matrix. In this work we 
have shown that one has a chance to identify the sign of 
through an oscillation parameter fitting performed 
in representation T). We found two local minima when 
analyzing global neutrino data in representation P, one 
corresponding to 6*^-^ < (the global minimum) and the 
other one efr^ > (the second minimum). A weak pref- 
erence for negative sin^ia is found. It is interesting to 
note that the global solar data point to > while 
all the other data < 0. We do not have interpreta- 
tion to this difference yet. One may investigate whether 
the matter (MSW) effects have impacts on the sign of 
Of^ during the propagation of solar neutrinos through 
the Sun and the Earth. This issue will be further stud- 
ied in a separate work. In addition, we also provide the 
information for the correlation of any two mixing angles. 
It is found that the correlations in representation A are 
less than those in 7), but the three mixing angles are not 
completely decoupled in representation A. 

In conclusion, owing to the strong correlations among 
the three mixing angles in the new parametrization, the 
advantages of doing the global neutrino oscillation anal- 
ysis using data from past, current, and near future neu- 
trino oscillation experiments shall become manifest. 
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Appendix A: Solar Neutrino Sector 

1. CI and Ga Rate Neutrino Data 

The rates for the Chlorine and Gallium experiments 
are expressed in terms of SNUs (1 Solar Neutrino Unit — 
one interaction per 1036 target atoms per second). The 
predicted rate for the Chlorine/ Gallium experiment is 
given by 

/>oo 

Rci/Ga = / dE,cl>,^{E,)P,,{E,)aci/Ga{E,), (Al) 

where E^, and Eth are respectively the neutrino energy 
and the threshold energy for neutrinos captured by chlo- 
rine or gallium, (j)^^ (E^) is the flux of electron neutrinos 
arriving at the detector, including the flux from all solar 
neutrino reactions, Pee{E^) is the survival probability of 
i^e — >■ i/e, (j{E^) is the cross section of electron neutrino 
with the target (either chlorine or gallium). The CI rate 
of 2.56 ± 0./6(statistical) ± 0.16(systematic) SNU from 
Ref. dl and the Ga rate of 66.1 ± 3.1 SNU from SAGE, 
Gallex, and GNO [i^l are used in this analysis. 



2. Borexino Rate Neutrino Data 
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The predicted rate calculation of ^Be solar neutrinos 
measured in Borexino experiment is calculated as 

Rbcc = i?L {SeiE,,E,)p, + S.AE.,E,){l-p,)} , 



where Rg 



(A2) 

is the expected rate for non-oscillated solar 
i^e that is 74 ± 4 counts/(day x 100 ton), Se,x{E^,Ee) 
describe the probability that the elastic scattering of a 
Ve.x of energy E^, with electrons produces a recoil electron 
of energy E^, and Pe is the survival probability of Ve — >■ 
Ve- The measured rate of 49 ± 3stat ± ^^vst counts/ (day 
• 100 ton) from Borexino experiment f48!| is adopted in 
this analysis. 



3. Super-Kamiokande Solar Neutrino Data 

The measured day/ni ght spectra of Super-Kamiokande 
(SK) phases I & II [ll, Il5| are employed in this analy- 
sis. To fit for the oscillation parameters, we follow the 
methods described in Ref. [l3| to build up the predicted 
rate calculation for SK solar day /night spectra, which is 
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given by 

RsK ^ I dE I dE^I{E^) I dEeR{Ee,E)x 

Jeo Je„ Je^ 

{Se{E,, E,)p, + S^,{E,,E,){1 - pe)} , (A3) 

where I{E^) is the spectrum of or hep neutrino, 
R{Ee,E) is the detector response function representing 
the probabihty that a recoil electron of energy E^ is re- 
constructed with energy E, and Se,x{Eu, E^) and Pe have 
the same meanings as given in previous paragraph. Both 
rates due to *B and hep interactions are taken into con- 
sideration in this work. 

4. SNO Solar Neutrino Data 

For SNO solar data, the fraction of the extracted CC 
(charged current) spectra as a function of one unoscil- 
lated SSM 5^ (Fig. 28 in Ref. ^) is used for the mixing 
parameter fitting. The theoretical calculation of the CC 
flux fraction for a set of oscillation parameter in order for 
comparison with SNO solar data is described as follows. 
The CC fiux fraction in the electron kinetic energy be- 
tween Ti and T2, denoted as f{Ti,T2), is calculated as 
the number of events observed in the data set contributed 
to CC interactions by the signal extraction with electron 
kinetic energies between Ti and T2 divided by the num- 
ber of all CC events that would be observed above the 
threshold kinetic energy, Tth- It is formulated as [sof 

/(Ti,r2) = 

/o°° /o°° 4>{E,)Pee{E,)^^{E,, T,)R{T,,T'^dE,dT,dTl 
r r /tI ^{E,)^^{El,T,)R{T,X)dE,dT,dT: 

(A4) 

where P^p is the survival probability of ^ Ve, (f'iP'u) in- 
cludes both fiuxes of '^B and hep solar neutrinos as a func- 
tion of neutrino energy, dcr/dTg is the differential cross 
section for CC interactions, and R{Te,Tp) is the energy 
resolution function, describing the probability of seeing 
an apparent kinetic energy Tg for a given true energy Te, 
whose expression can be found in Ref. [l^ . 



5. Global Solar Neutrino Data 

To fit for oscillation parameters using global solar data, 
the combined for the observables and predictions is 
given by 

Xso/ar 
N 

J2 (O. - Or)[4-(tot)]-i(0, - Of P), (A5) 

where i, j are indices to sum over all the observables. 
The quantities Oi 

Oi = ©^(ATO^tan^e,...) 



is the theoretical prediction for that observable, and Of p 
represents a series of observables from a number of so- 
lar neutrino experiments that include rate measurements 
from CI and Ga (e.g. Homestake, Gallex, GNO, and 
SAGE) experiments, rate of ''Be solar neutrino measure- 
ment from Borexino experiment, a number of spectral 
shape measurements from SK and SNO as mentioned 
above. The total error matrix, (T^(toi), is a sum of contri- 
butions from the rate and spectral measurements, which 
is given as 



a'^{tot) = 
cr^ {exp) + (j\ 



(y\{Bx) + aliSNO) + aUSK), (A6) 



where a'^{exp) is a diagonal matrix containing the statis- 
tical and systematic errors from the rate measurements 
and the statistical errors from the spectral measurements; 
ct|j is the rate error matrix, handling the correlations be- 
tween the rate CL and Ga experiments using the proce- 
dure described in Ref. [60*1; cr^(i?x) is treated as indepen- 
dent from cr^ in our analysis. Furthermore, the spectral 
correlation matrix a'c;{SNO) for SNO measurements is 
assumed to be uncorrelated with (Tg{SK) for SK mea- 
surements, which are uncorrelated with the CL, Ga, and 
Borexino measurements. One can refer to Refs. |60l - [63j 
for detailed discussion on the covariance error matrix, 
a^(tot)^^. By minimizing xioiar with three- fiavor neu- 
trino oscillation scheme while fixing Am^j = 2.4 x 10"'^ 
eV^ and sin^ 623 =0.5, the values of Amji, tan^ 9i2, and 
sin^ 9i3 were found to reproduce the results of global so- 
lar data as reported in Ref. [l^ when the conventional 
mixing matrix parametrization is used. 



Appendix B: Reactor Neutrino Sector 

1. KamLAND Reactor Data 

For KamLAND reactor data, the prompt energy spec- 
trum of Ve, events of energy (as shown in Fig. (1) of 



Ref. 
Ref. 



16 1) is used. We follow the strategy described in 
64l | for this analysis. Since the distance from the 
reactor source to the KamLAND detector is ^^180 Km, 
the treatment of vacuum oscillation is considered in this 
analysis. With a little modification, the number of ex- 
pected events per unit of the prompt position energy is 
given by 



NiT:)=^At-e{T:)- J dE^^^ 



(Bl) 



The values $ = 5.98 x lO^^ protons and At = 2135 days 
are the total number of target protons and livetime, re- 
spectively. e(Tg) is the detection efficiency as a function 
of measured prompt energy, . d(j)/dE^ is the time- 
averaged differential neutrino flux at the KamLAND de- 
tector, where the flssion flux as a function of distance can 
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be found at the website (65| and the relative fission yields, 
235u.238u.239py.24ip^^ is taken from [l6|, and flux from 

Korean reactors is not taken into account since only ~3% 
of the total flux contributes to KamLAND signal. In the 
presence of oscillation, each jth reactor term in d(f>/dEi, 
must be multiplied by the corresponding neutrino sur- 
vival probability Pee{Eu, Lj) for neutrinos of energy E"^ 
and taveling distance of Lj between each reactor source 
and the KamLAND. i?(Te,Tg) is the energy resolution 
function with Gaussian width equal to QA%yjTjWeV , 
which has the same meaning as described in the solar 
sector. da(Ej^,Te)/dTe is the inverse beta decay cross 
section (63.l66j. 

To fit for the oscillation parameters, the function 
is of the Gaussian form [g^I that includes both the total 
number of events and the spectral shape: 



were used. These results consist of seven positron energy 
bins for each of the two reactors for a total of 14 bins. The 



covariant matrix ■ defined in equation 54 in Ref. jl7| 



is applied to the function. This matrix accounts for 
any correlations between the energy bins of the two re- 
actors. We then minimize the function with respect 
to the neutrino oscillation parametrization and we can 



reproduce the CHOOZ results as presented in Ref. [17 1 
when the conventional mixing matrix parametization is 
applied. 



Appendix C: Long-Baseline Accelerator Neutrino 
Sector 

1. MINOS 



Xkl 



Ajtheo 



jvrobs 
'-^tot 



E 



N. 



theo 



obs 



(B2) 



where N°^^ — 1780 is the total number of observed events 
after substracting all background events, Nf^^" is the pre- 
dicted total events; iV*'''=° and N"'"' are the numbers of 
observed events and predicted events respectively in zth 
energy bin, the total error, cr^^**^, is the quadrature sum 
of the statistical and systematic uncertainties, and (t[°*'^ 
is the error for ith bin by adding statistical and system- 
atic uncertainties in quadrature. Minimizing the Xkl 
while fixing Amgj = 2.4 x 10^^ eV^ and sin^ 623 = 0.5, 
the values for Am2i, tan^6'i2, and sin^ 6113 were found 
to reproduce the KamLAND results in Ref. [l6| if the 
conventional mixing matrix parametrization is applied. 



CHOOZ Reactor i^e Data 



The analysis for MINOS used results from the i/^ disap- 
pearance channel and the i^e appearance channel [20| . 
Results from both channels consisted of the reconstructed 
neutrino energy as seen in the far detector. The distance 
of 735 km between the near and far detectors was con- 
sidered to be too small for matter effects to have any 
significant contribution to the results. Vacuum oscilla- 
tion is thus applied in this analysis. 

The predicted number of neutrinos at the far detector 
is calculated as: 



N" 



NC 



(CI) 



where P^^ is the survival probability for j/^, P^r is the 
probability for Vr appearance, is the Vt detection effi- 
ciency, N'^'^ is the NC (neutral current) spectrum, / is 
the normalization factor, and g is the NC scaling factor. 
Like the data, the data consists of events binned 
by reconstructed energy over a range of 1 to 8 GeV. The 
predicted neutrino spectrum of i^e appearance channel is 
calculated using: 



The CHOOZ results [T3| reported the observed 
positron energy from the neutrino interactions. For this 
experiment, the distance between the detector and the 
neutrino source is relatively small (1 km) and thus the 
neutrino propagation through vacuum can be used to cal- 
culate the anticipated positron spectrum. Using the pro- 
cedures reported in Ref. [13] , the expected positron yield 
for the fc-th reactor and the j-th energy spectrum bin is 
parametrized as 



-X{E,) 
U - 1, 



P{E,,Lk,9,Sml2), 
..,7, fc = l,2), (B3) 



where X{Ej) is the distance-independent positron yield 
for no presence of neutrino oscillation, ife is the reactor- 
detector distance, and P{Ej,Lk, 9, 6171^2) is the oscilla- 
tion probability averaged over the energy bin and the 
core sizes of the detector and the reactor. For the exper- 
imental data, the results reported in Table 8 in Ref. [l3| 



N° 



gN 



NC 



(C2) 



where P^e is the probability for i^e appearance, and / 
{g ) has the same meaning as / (g) but with a different 
value. The value for iV"° is taken from [19]. 

The function is of the Poisson form and includes 
both the spectral shape and the total number of neutrino 
events, which is given by: 



E 



ln(iV-7iV°''/))] 
iVf^(l + ln(iVr7A))] 



(C3) 



where iV°„f and 7V°j'/ 



are the predicted and detected to- 
tal number of neutrinos respectively; TV/"*^ and iV°'"' are 
the predicted and detected number of neutrinos in the 
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i energy bin. The systematics for the energy scale, 
normalization factor and scale factor for the NC events 
are represented by with uncertainty of aj. Minimiz- 
ing the foi' the j/^ disappearance channel while fixing 
Amii = 7.67 X 10"^ eV^, tan^ 612 = 0.427, sin^ 6*13 = 0.0 
and assuming a normal mass hierarchy, the values for 
Am§2 and sin^^23 where found to reproduce the MINOS 
results in Ref. [l9[ when the mixing matrix parametriza- 
tion A is applied. 



2. K2K 

The results from K2K [l^ used in this analysis were 
the reconstructed energy spectrum for one-ring fi- 
like sample. Like MINOS, the K2K's 250-km oscihation 
length is considered to be too short to manifest matter 
enhanced oscillation effects. Vacuum oscillation can thus 
be applied to this analysis. 

The predicted no-oscillation neutrino spectrum rep- 
resents the true neutrino energy spectrum at the near 
detector. To extract the predicted far detector oscilla- 
tion spectrum from the no-oscillation spectrum, the neu- 
trino interaction cross sections and detection efficiencies 
must be applied. Furthermore, both CC and NC events 
are present in the neutrino spectrum and hence are ac- 
counted for separately. The expected number of neutrino 
events for oscillation neutrinos is given by 

^osc ^ ^no osc [g^^p^^^^^ + ,^^^^^] (C4) 

where CTcc and cr„c are the cross sections for the CC and 
NC interactions respectively, Ccc and e„c are the detec- 
tion efficiencies for the CC and NC interactions, and Pec 
is the survival probability for v^. The energy response 
function is applied to obtain the probability of seeing the 
measured energy for a given true energy, same meaning 
as employed in the solar analyses. The form of the energy 
response function is employed from Ref. |67|: 



Appendix D: Atmospheric Neutrino Sector 

1. Super-Kamiokande Atmospheric Data 

The Super-Kamiokande (SK) detector is located deep 
under the peak of Mt Ikenoyama, where the 1,200-m rock 
overburden can reduce the flux of cosmic rays reaching 
the detector down to about 3 Hz (see e.g. Ref. 68]). 
Atmospheric neutrinos penetrating the Earth interact 
with the nucleus or nucleons in the SK water tank or in 
the surrounding rock, giving rise to partially contained 
(PC), fully contained (FC), or upward-going muon (UP/i) 
events. Two and three-flavour neutrino oscillation anal- 
yses [21I, m, [6^ have been performed using the accumu- 
lated data from SK-I, II and III phases. 

We try to reproduce the results of the 3-flavour anal- 
ysis published in Ref. [2^ to include in our global anal- 
ysis. Therein, the 1489 live-days of PC and FC, and 
the 1646 days UP^ data collected during the SK-I pe- 
riod (1996-2001) are used. Based on the types of the 
out-going leptons, their energy deposited in the SK de- 
tector, and their zenith angle (— 1 < cos^zcnith < 1 for 
PC and FC; —1 < cos^zonith < for UP/i), the selected 
events are divided in 370 bins. In our work, however, 
lack of information forced us to follow the approach of 
Ref. [67| . where only 55 energy and zenith angle bins of 
SK atmospheric data are used. 

Here we briefly describe our analysis. First, the NU- 
ANCE package [53] is used to simulate atmospheric neu- 
trino events assuming no oscillation effects. For achieving 
good statistics we have run the simulations for 200-year 
operation time. Neutrino oscillation effects in the atmo- 
sphere and inside the Earth are then incorporated by the 
"weighting" factors [68] as: 

We = P{Ve Ve) + ^ '^e) , for , 

Wf, = P(i/p -> v^) + ^P{ve i^/.) , for v^, ■ (Dl) 

YIJ. 



N 



obs 



' N? 



* e 



-i((£;„6,-£;,+o.05)/o-)=^ 



(C5) 



where Eobs and Ei are the measured and true neutrino 
energies respectively, and cr is defined as p^: 



cr = 0.2Ei{l - e 



[Q.2-Ei)/Q.B 



The K2K analysis uses the same function as that 
used in the MINOS analysis. By minimizing the for 
the while fixing Am|i = 7.67 x 10"^ cV^, tan^ 6*12 = 
0.427, sin^ ^13 = 0.0, and assuming a normal mass hi- 
erarchy, the values for Arrigj and sin^ 623 were found to 
reproduce the K2K results in Ref. [l8| when the conven- 
tional mixing parametrization is applied. 



For the (energy and angular-dependent) incident atmo- 
spheric neutrino fluxes 4>e,^n we also adopt the Honda 
three-dimensional calculation [70]. We take into account 
the effects when neutrinos oscillating in matter following 
the prescription of Ref. |58| , and approximate the Earth 
as a four-layer division as Ref. [23 1 did. Each layer has 
a constant density (inner core: R < 1221 km, p = 13.0 
g/cm^; outer core: 1221 < R < 3480 km, p = 11.3 
g/cm^; mantle: 3480 < R < 5701 km, p = 5.0 g/cm^; 
crust: 5701 < R < 6371 km, p = 3.3 g/cm^). We ap- 
ply similar criteria and cuts on the kinematics of the 
simulated events so that we achieve the same selection 
efficiency as Ref. [2^ does. 

Furthermore, one can include the systematics by using 
the " pull technique" [i^, |7l| . Due to correlated system- 
atic uncertainties, the event rate in the n-th bin predicted 
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by the MC simulation, i?*^*^, is shifted by an amount as 



R 



MC 



MC 



' 11 ^ 



Here S^^ is the Icr error associated to the fc-th source 
of systematics, and ^^'s are a set of univariate Gaussian 
random variables. Ref. (67j gives a summary and detailed 
discussion of the 11 systematics sources, which we also 
use in our work. The function thus becomes 



XSK — 



Ma-} 




k,h=l 



(D2) 



regions for 100 year operation time, assuming no oscil- 
lation. Oscillation effects are then added with the help 
of the "weighting" factor. For the incident atmospheric 
neutrino fluxes, we adopt the Bartol three-dimensional 
calculation [t^]. Due to lack of a detector simulation 
package such as the SNOMAN, we apply simple cuts on 
the kinematics of simulated events, separately for differ- 
ent event types (f^-induced through-going muons, u^^ wa- 
ter interactions, and i^f^ and Ug internal interactions). We 
require (i) the impacat parameter b < 830 cm; (ii) the 
muon energy when entering the detector > 1 GeV, 
and adjust the "trigger efhciency" so that the yearly 
event rate match that given in Ref. 22]. The trigger 
efficiency we defined summarises all other instrumental 
cuts which we are not able to apply. In addition, the 
systematics are taken into account using the generalised 
" pull technique" . Following Ref. [l^, [t^] , the likelihood 
function used in our analysis is 



for the observed event rate R'^ in each bin. Here cr**°* is 
the statistical error of the n-th bin, and the inverse 
of the correlation matrix, the value of which can also be 
found in Ref. [13 • 

We minimize the Xsk value with respect to the oscilla- 
tion parameters by solving — 0, which is equivalent 
to solving a set of linear equations: 



11 

E 



fc=l Ln=l 



55 / dMC \ 2 



55 

E 

n=l 
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(D3) 



Jtotal 



= - In/: 



2E 



N: 



data 
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data 



^data^ 



(D4) 



where TV^'^ata ^j^g measured event number in zenith an- 
gle bin j, and N^^f denotes the expected event number 
in bin i with the systematics a equal to zero. We include 
8 systematics error sources, and calculate their values 
which minimise the likelihood function 



Ea 



r 



data 



5" 



(D5) 



Under normal hierarchy and A parametrization, our 
best-fit values for (Amgj, sin^023, sin^^ia) are (2.4 x 
10^3 0.45, 0.0), and the 90% confidence levels 



4.61) are 1.5 x 10"'^ eV^ < 



10 ev , u.ou ^ sin (723 <k u.uy, auu sm C13 
which is consistent with those published in Ref. [2 



Am|2 < 3.3 X 

3 eV^ 0.35 < SIT? 62?. < 0.59, and sin^ 6*13 < 0.17, 



Here the coefficients /3i are available in Ref. [7^ , and the 
new error matrix is 



(D6) 



2. SNO Atmospheric Data 

Thanks to its deep location and flat overburden, the 
SNO detector can observe atmospheric neutrinos over a 
wide range of zenith angle, via their charged-current in- 
teractions in the surrounding rock. Angular distribution 
of through-going muons having —1 < cos6'zonith < 0.4 
can be used to infer neutrino oscillation parameters and 
incident neutrino flux, while data above this cutoff pro- 
vide access to the study of the cosmic-ray muon flux. 

Below we describe briefly our approach to including the 
latest SNO atmospheric neutrino data from Ref. [13] . De- 
tails of the analyses can be found in T. Sonle y's PhD the- 
sis j^]. We first use the NUANCE package |57| to simu- 
late atmospheric neutrino induced through-going muons 
and muons generated in SNO detector's D2O and H2O 



with a 2 the diagonal error matrix whose entries repre- 
sent the size of the systematic errors. 

We first compare the results of our two-flavour anal- 
ysis with those in Ref. [12, [T^I. Normalisation of the 
Bartol three-dimensional atmospheric neutrino flux $0 is 
determined simultaneously with the neutrino oscillation 
parameters. Our best flt values for {Am^2 7 sin'^ 2023 , $0) 
are (2. Ox 10^'^ , 1.0 , 1.30) in the normal hierarchy scheme. 
These are close to those obtained in Ref. [t^, while still 
within Icr of those from Ref. H^]. To save CPU time, we 
then do the three-flavour analysis by keeping $0 flxed at 
1.30. Matter effects induced when neutrinos pass through 
different Earth la yer s are considered following the pre- 
scription of Ref. [5^. As expected, SNO atmospheric 
neutrino data are not sensitvie to 613. Since, unlike 
the Super-Kamiokande experiment, SNO only observes 
muons plus a few possible i/g-induced internal events. 



Appendix E: 



9^,5^) Solutions in Terms of {9fj,5fp) 



Following the procedure presented in Ref. j9|, the 
transformation of {Oj^,S^p) from representation 23 to ^ 
is again briefly summarized here by only listing the solu- 
tions of the nine parameters. 

To solve for {d:^,S'^) in representation A when the 



,S^p) are known in representation 2?, begin with 
U = RWR{ef^, S^^ = D"^ ■ RWRief^, SZ) ■ D"" , (El) 



where 



cos (7i ^ 



'13 



13 



(E8) 



-12 



The remaining parameters thus can be determined as fol- 
lows. 



sin' $m 



(sin 



L3 



risin<5g)^ 



(sin - dsmSfp)^ + (d cos - cos $^3) 



2 ' 



(E9) 



(E2) 



Through the nine real parts and the nine imaginary parts 
of Equation ([E]) , the solutions to the nine parameters are 
listed as follows: 



sin' $;^3 



*L2 = 0, 
*Li + *m = 0. 



(a sin (55,) 2 



(a sinS'gy + (a cos 62, - by 



where a — S23'Si2'Si3 and b — (^23(^13 
sin^ $^ 



(a sin (55,) 2 



(a' sin (5?,) 2 -|- (a' cos 55, — b') 



(E3) 



(E4) 



(E5) 



where a 



„D „D „D I. _ -D -D J 
C23'^12Cl3; " — ■S23'5l3' 'i'^g 



LI 



*^fl2- With these phases, the three mixing angles can be 
extracted in representation A: 



s?3cf2COS$^2 



tan6'5^ = 



ciscfa cos 



L3 



V V T> 
*23'^12'5l3 



4cos(<i>^3~(55J 



COS! 



cf3Sl2cf3 C0S($ 



12 



(5^) 



S?3S?3COS$4; 



15 



:;f2cf3COs($^i + $^i) 



(E6) 



(E7) 



where 



'-23*12 



<ffl 



. *23''12*13 

and ^i. 



can be 



Therefore 

determined. Finally, the CP-violating phase in repre- 
sentation A, 6^, can be resolved using those conditions 
associated with sin (5^: 



sin (5^ 



L2 



o*-^ 

*23''12*13 



f2sf3sin($^3 + 'fm) 



/-i-^ /-i-^ 

''23'-12*13 



<i>^2) 



*23'^12'^13 



23^i2^f3 sin(?? 



'^23*12''13 



23^12^f3 Shl($fi 



"-23* 13 



sin$ 



Ll 



(ElO) 



'13 



where r/ = $-^3 -I- $7^2- Consistency in the value of s in 6-^ 
calculated from the five different expressions in Eq. (jElOp 
serves as a means to check whether the values of the nine 
parameters are correct. 
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